! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 
! The ODE Function of Chemical Model File
! 
! Generated by KPP-2.2.3 symbolic chemistry Kinetics PreProcessor
!       (http://www.cs.vt.edu/~asandu/Software/KPP)
! KPP is distributed under GPL, the general public licence
!       (http://www.gnu.org/copyleft/gpl.html)
! (C) 1995-1997, V. Damian & A. Sandu, CGRER, Univ. Iowa
! (C) 1997-2005, A. Sandu, Michigan Tech, Virginia Tech
!     With important contributions from:
!        M. Damian, Villanova University, USA
!        R. Sander, Max-Planck Institute for Chemistry, Mainz, Germany
! 
! File                 : aromatics_kpp_Function.f90
! Time                 : Thu Jan  7 01:53:05 2021
! Working directory    : /n/home08/kbates/Aromatics/GC_new3
! Equation file        : aromatics_kpp.kpp
! Output root filename : aromatics_kpp
! 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~



MODULE aromatics_kpp_Function

  USE aromatics_kpp_Parameters
  IMPLICIT NONE

! A - Rate for each equation
  REAL(kind=dp) :: A(NREACT)

CONTAINS


! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 
! Fun - time derivatives of variables - Agregate form
!   Arguments :
!      V         - Concentrations of variable species (local)
!      F         - Concentrations of fixed species (local)
!      RCT       - Rate constants (local)
!      Vdot      - Time derivative of variable species concentrations
! 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

SUBROUTINE Fun ( V, F, RCT, Vdot )

! V - Concentrations of variable species (local)
  REAL(kind=dp) :: V(NVAR)
! F - Concentrations of fixed species (local)
  REAL(kind=dp) :: F(NFIX)
! RCT - Rate constants (local)
  REAL(kind=dp) :: RCT(NREACT)
! Vdot - Time derivative of variable species concentrations
  REAL(kind=dp) :: Vdot(NVAR)


! Computation of equation rates
  A(1) = RCT(1)*V(69)
  A(2) = RCT(2)*V(69)*V(169)
  A(3) = RCT(3)*V(69)*V(164)
  A(4) = RCT(4)*V(69)*V(165)
  A(5) = RCT(5)*V(69)*V(165)
  A(6) = RCT(6)*V(34)
  A(7) = RCT(7)*V(164)*V(169)
  A(8) = RCT(8)*V(165)*V(169)
  A(9) = RCT(9)*V(164)*V(164)
  A(10) = RCT(10)*V(164)*V(166)
  A(11) = RCT(11)*V(165)*V(166)
  A(12) = RCT(12)*V(165)*V(166)
  A(13) = 5.35e+07*V(34)
  A(14) = RCT(14)*V(169)*V(175)
  A(15) = RCT(15)*V(64)*V(175)
  A(16) = RCT(16)*V(73)*V(175)
  A(17) = RCT(17)*V(169)*V(172)
  A(18) = RCT(18)*V(172)*V(175)
  A(19) = RCT(19)*V(172)*V(172)
  A(20) = RCT(20)*V(164)*V(175)
  A(21) = RCT(21)*V(165)*V(175)
  A(22) = 2e-11*V(166)*V(175)
  A(23) = RCT(23)*V(164)*V(172)
  A(24) = RCT(24)*V(165)*V(172)
  A(25) = RCT(25)*V(58)*V(175)
  A(26) = 4e-12*V(166)*V(172)
  A(27) = RCT(27)*V(50)*V(175)
  A(28) = RCT(28)*V(117)*V(175)
  A(29) = RCT(29)*V(169)
  A(30) = RCT(30)*V(169)
  A(31) = RCT(31)*V(73)
  A(32) = RCT(32)*V(165)
  A(33) = RCT(33)*V(166)
  A(34) = RCT(34)*V(166)
  A(35) = RCT(35)*V(50)
  A(36) = RCT(36)*V(117)
  A(37) = RCT(37)*V(166)*V(166)
  A(38) = RCT(38)*V(38)
  A(39) = RCT(39)*V(58)
  A(40) = RCT(40)*V(38)
  A(41) = RCT(41)*V(58)
  A(42) = RCT(42)*V(42)*V(175)
  A(43) = RCT(43)*V(164)*V(170)
  A(44) = RCT(44)*V(170)*V(172)
  A(45) = RCT(45)*V(170)*V(170)
  A(46) = RCT(46)*V(170)*V(170)
  A(47) = RCT(47)*V(165)*V(170)
  A(48) = RCT(48)*V(52)
  A(49) = RCT(49)*V(57)*V(175)
  A(50) = RCT(50)*V(57)*V(175)
  A(51) = RCT(51)*V(149)*V(175)
  A(52) = RCT(52)*V(149)*V(166)
  A(53) = RCT(53)*V(122)*V(175)
  A(54) = RCT(54)*V(173)*V(175)
  A(55) = RCT(55)*V(57)
  A(56) = RCT(56)*V(149)
  A(57) = RCT(57)*V(149)
  A(58) = RCT(58)*V(52)
  A(59) = RCT(59)*V(52)
  A(60) = RCT(60)*V(35)*V(175)
  A(61) = RCT(61)*V(164)*V(171)
  A(62) = RCT(62)*V(171)*V(172)
  A(63) = RCT(63)*V(171)*V(171)
  A(64) = RCT(64)*V(171)*V(171)
  A(65) = RCT(65)*V(170)*V(171)
  A(66) = RCT(66)*V(170)*V(171)
  A(67) = RCT(67)*V(165)*V(171)
  A(68) = RCT(68)*V(56)
  A(69) = RCT(69)*V(61)*V(175)
  A(70) = RCT(70)*V(61)*V(175)
  A(71) = RCT(71)*V(126)*V(175)
  A(72) = RCT(72)*V(126)*V(166)
  A(73) = RCT(73)*V(99)*V(175)
  A(74) = RCT(74)*V(43)*V(175)
  A(75) = RCT(75)*V(61)
  A(76) = RCT(76)*V(126)
  A(77) = RCT(77)*V(126)
  A(78) = RCT(78)*V(56)
  A(79) = RCT(79)*V(56)
  A(80) = RCT(80)*V(128)*V(171)
  A(81) = RCT(81)*V(144)*V(171)
  A(82) = RCT(82)*V(147)*V(171)
  A(83) = RCT(83)*V(158)*V(171)
  A(84) = RCT(84)*V(162)*V(171)
  A(85) = RCT(85)*V(142)*V(171)
  A(86) = RCT(86)*V(148)*V(171)
  A(87) = RCT(87)*V(137)*V(171)
  A(88) = RCT(88)*V(145)*V(171)
  A(89) = RCT(89)*V(161)*V(171)
  A(90) = RCT(90)*V(161)*V(171)
  A(91) = RCT(91)*V(167)*V(171)
  A(92) = RCT(92)*V(167)*V(171)
  A(93) = RCT(93)*V(78)*V(175)
  A(94) = RCT(94)*V(78)*V(175)
  A(95) = RCT(95)*V(91)*V(175)
  A(96) = RCT(96)*V(159)*V(175)
  A(97) = RCT(97)*V(159)*V(166)
  A(98) = RCT(98)*V(165)*V(167)
  A(99) = RCT(99)*V(67)
  A(100) = RCT(100)*V(164)*V(167)
  A(101) = RCT(101)*V(47)*V(175)
  A(102) = RCT(102)*V(147)*V(164)
  A(103) = RCT(103)*V(36)*V(175)
  A(104) = RCT(104)*V(36)*V(175)
  A(105) = RCT(105)*V(128)*V(164)
  A(106) = RCT(106)*V(144)*V(164)
  A(107) = RCT(107)*V(45)*V(175)
  A(108) = RCT(108)*V(158)*V(164)
  A(109) = RCT(109)*V(158)*V(164)
  A(110) = RCT(110)*V(162)*V(164)
  A(111) = RCT(111)*V(142)*V(164)
  A(112) = RCT(112)*V(148)*V(164)
  A(113) = RCT(113)*V(154)*V(164)
  A(114) = RCT(114)*V(154)
  A(115) = RCT(115)*V(141)*V(164)
  A(116) = RCT(116)*V(152)*V(164)
  A(117) = RCT(117)*V(152)
  A(118) = RCT(118)*V(132)*V(164)
  A(119) = RCT(119)*V(137)*V(164)
  A(120) = RCT(120)*V(157)*V(164)
  A(121) = RCT(121)*V(157)*V(166)
  A(122) = RCT(122)*V(150)*V(166)
  A(123) = RCT(123)*V(135)*V(166)
  A(124) = RCT(124)*V(135)*V(164)
  A(125) = RCT(125)*V(135)*V(165)
  A(126) = RCT(126)*V(135)*V(172)
  A(127) = RCT(127)*V(123)*V(166)
  A(128) = RCT(128)*V(123)*V(164)
  A(129) = RCT(129)*V(123)*V(172)
  A(130) = RCT(130)*V(123)*V(170)
  A(131) = RCT(131)*V(150)*V(169)
  A(132) = RCT(132)*V(150)*V(175)
  A(133) = RCT(133)*V(120)*V(164)
  A(134) = RCT(134)*V(120)*V(172)
  A(135) = RCT(135)*V(120)*V(170)
  A(136) = RCT(136)*V(157)*V(157)
  A(137) = RCT(137)*V(145)*V(164)
  A(138) = RCT(138)*V(45)*V(166)
  A(139) = RCT(139)*V(174)*V(175)
  A(140) = RCT(140)*V(116)*V(175)
  A(141) = RCT(141)*V(163)*V(175)
  A(142) = RCT(142)*V(161)*V(165)
  A(143) = RCT(143)*V(40)
  A(144) = RCT(144)*V(55)*V(165)
  A(145) = RCT(145)*V(146)*V(165)
  A(146) = RCT(146)*V(110)
  A(147) = RCT(147)*V(100)
  A(148) = RCT(148)*V(161)*V(164)
  A(149) = RCT(149)*V(146)*V(164)
  A(150) = RCT(150)*V(163)*V(166)
  A(151) = RCT(151)*V(108)*V(175)
  A(152) = RCT(152)*V(128)*V(170)
  A(153) = RCT(153)*V(144)*V(170)
  A(154) = RCT(154)*V(158)*V(172)
  A(155) = RCT(155)*V(162)*V(172)
  A(156) = RCT(156)*V(142)*V(172)
  A(157) = RCT(157)*V(148)*V(172)
  A(158) = RCT(158)*V(154)*V(172)
  A(159) = RCT(159)*V(141)*V(172)
  A(160) = RCT(160)*V(152)*V(172)
  A(161) = RCT(161)*V(132)*V(172)
  A(162) = RCT(162)*V(137)*V(172)
  A(163) = RCT(163)*V(157)*V(172)
  A(164) = RCT(164)*V(145)*V(172)
  A(165) = RCT(165)*V(168)*V(175)
  A(166) = RCT(166)*V(147)*V(170)
  A(167) = RCT(167)*V(166)*V(168)
  A(168) = RCT(168)*V(158)*V(170)
  A(169) = RCT(169)*V(162)*V(170)
  A(170) = RCT(170)*V(142)*V(170)
  A(171) = RCT(171)*V(148)*V(170)
  A(172) = RCT(172)*V(154)*V(170)
  A(173) = RCT(173)*V(154)*V(154)
  A(174) = RCT(174)*V(143)*V(170)
  A(175) = RCT(175)*V(141)*V(170)
  A(176) = RCT(176)*V(152)*V(170)
  A(177) = RCT(177)*V(132)*V(170)
  A(178) = RCT(178)*V(137)*V(170)
  A(179) = RCT(179)*V(157)*V(170)
  A(180) = RCT(180)*V(145)*V(170)
  A(181) = RCT(181)*V(72)*V(175)
  A(182) = RCT(182)*V(106)*V(175)
  A(183) = RCT(183)*V(147)*V(147)
  A(184) = RCT(184)*V(147)*V(147)
  A(185) = RCT(185)*V(147)*V(172)
  A(186) = RCT(186)*V(128)*V(172)
  A(187) = RCT(187)*V(144)*V(172)
  A(188) = RCT(188)*V(167)*V(172)
  A(189) = RCT(189)*V(161)*V(172)
  A(190) = RCT(190)*V(146)*V(172)
  A(191) = RCT(191)*V(115)*V(175)
  A(192) = RCT(192)*V(115)*V(169)
  A(193) = RCT(193)*V(100)*V(175)
  A(194) = RCT(194)*V(110)*V(175)
  A(195) = RCT(195)*V(100)*V(169)
  A(196) = RCT(196)*V(4)*V(175)
  A(197) = RCT(197)*V(138)*V(175)
  A(198) = RCT(198)*V(138)*V(175)
  A(199) = RCT(199)*V(115)*V(166)
  A(200) = RCT(200)*V(125)*V(175)
  A(201) = RCT(201)*V(151)*V(175)
  A(202) = RCT(202)*V(125)*V(166)
  A(203) = RCT(203)*V(151)*V(166)
  A(204) = RCT(204)*V(114)*V(175)
  A(205) = RCT(205)*V(160)*V(175)
  A(206) = RCT(206)*V(156)*V(175)
  A(207) = RCT(207)*V(155)*V(175)
  A(208) = RCT(208)*V(155)*V(175)
  A(209) = RCT(209)*V(128)*V(167)
  A(210) = RCT(210)*V(144)*V(167)
  A(211) = RCT(211)*V(128)*V(167)
  A(212) = RCT(212)*V(144)*V(167)
  A(213) = RCT(213)*V(114)*V(169)
  A(214) = RCT(214)*V(160)*V(169)
  A(215) = RCT(215)*V(156)*V(169)
  A(216) = RCT(216)*V(114)*V(166)
  A(217) = RCT(217)*V(156)*V(166)
  A(218) = RCT(218)*V(156)*V(166)
  A(219) = RCT(219)*V(161)*V(170)
  A(220) = RCT(220)*V(146)*V(170)
  A(221) = RCT(221)*V(161)*V(170)
  A(222) = RCT(222)*V(146)*V(170)
  A(223) = RCT(223)*V(98)*V(175)
  A(224) = RCT(224)*V(98)*V(175)
  A(225) = RCT(225)*V(71)*V(175)
  A(226) = RCT(226)*V(63)*V(175)
  A(227) = RCT(227)*V(59)*V(175)
  A(228) = RCT(228)*V(60)*V(175)
  A(229) = RCT(229)*V(75)*V(175)
  A(230) = RCT(230)*V(68)*V(175)
  A(231) = RCT(231)*V(70)*V(175)
  A(232) = RCT(232)*V(88)*V(175)
  A(233) = RCT(233)*V(81)*V(175)
  A(234) = RCT(234)*V(89)*V(175)
  A(235) = RCT(235)*V(88)*V(175)
  A(236) = RCT(236)*V(81)*V(175)
  A(237) = RCT(237)*V(89)*V(175)
  A(238) = RCT(238)*V(6)*V(175)
  A(239) = RCT(239)*V(48)*V(175)
  A(240) = RCT(240)*V(46)*V(175)
  A(241) = RCT(241)*V(41)*V(175)
  A(242) = RCT(242)*V(121)*V(172)
  A(243) = RCT(243)*V(121)*V(164)
  A(244) = RCT(244)*V(84)*V(175)
  A(245) = RCT(245)*V(83)*V(175)
  A(246) = RCT(246)*V(80)*V(175)
  A(247) = RCT(247)*V(80)*V(175)
  A(248) = RCT(248)*V(92)*V(175)
  A(249) = RCT(249)*V(92)*V(175)
  A(250) = RCT(250)*V(136)*V(167)
  A(251) = RCT(251)*V(136)*V(167)
  A(252) = RCT(252)*V(136)*V(170)
  A(253) = RCT(253)*V(136)*V(136)
  A(254) = RCT(254)*V(136)*V(172)
  A(255) = RCT(255)*V(136)*V(164)
  A(256) = RCT(256)*V(136)*V(164)
  A(257) = RCT(257)*V(49)*V(175)
  A(258) = RCT(258)*V(47)*V(166)
  A(259) = RCT(259)*V(153)*V(169)
  A(260) = RCT(260)*V(167)*V(167)
  A(261) = RCT(261)*V(167)*V(170)
  A(262) = RCT(262)*V(167)*V(170)
  A(263) = RCT(263)*V(158)*V(167)
  A(264) = RCT(264)*V(142)*V(167)
  A(265) = RCT(265)*V(148)*V(167)
  A(266) = RCT(266)*V(154)*V(167)
  A(267) = RCT(267)*V(143)*V(167)
  A(268) = RCT(268)*V(141)*V(167)
  A(269) = RCT(269)*V(152)*V(167)
  A(270) = RCT(270)*V(137)*V(167)
  A(271) = RCT(271)*V(162)*V(167)
  A(272) = RCT(272)*V(132)*V(167)
  A(273) = RCT(273)*V(157)*V(167)
  A(274) = RCT(274)*V(145)*V(167)
  A(275) = RCT(275)*V(158)*V(167)
  A(276) = RCT(276)*V(142)*V(167)
  A(277) = RCT(277)*V(148)*V(167)
  A(278) = RCT(278)*V(154)*V(167)
  A(279) = RCT(279)*V(143)*V(167)
  A(280) = RCT(280)*V(141)*V(167)
  A(281) = RCT(281)*V(152)*V(167)
  A(282) = RCT(282)*V(162)*V(167)
  A(283) = RCT(283)*V(132)*V(167)
  A(284) = RCT(284)*V(157)*V(167)
  A(285) = RCT(285)*V(145)*V(167)
  A(286) = RCT(286)*V(137)*V(167)
  A(287) = RCT(287)*V(147)*V(167)
  A(288) = RCT(288)*V(147)*V(167)
  A(289) = RCT(289)*V(161)*V(167)
  A(290) = RCT(290)*V(146)*V(167)
  A(291) = RCT(291)*V(112)*V(169)
  A(292) = RCT(292)*V(139)*V(169)
  A(293) = RCT(293)*V(153)*V(175)
  A(294) = RCT(294)*V(143)*V(164)
  A(295) = RCT(295)*V(143)*V(164)
  A(296) = RCT(296)*V(143)*V(172)
  A(297) = RCT(297)*V(139)*V(175)
  A(298) = RCT(298)*V(112)*V(175)
  A(299) = RCT(299)*V(130)*V(164)
  A(300) = RCT(300)*V(130)*V(172)
  A(301) = RCT(301)*V(134)*V(164)
  A(302) = RCT(302)*V(134)*V(172)
  A(303) = RCT(303)*V(90)*V(175)
  A(304) = RCT(304)*V(118)*V(175)
  A(305) = RCT(305)*V(131)*V(175)
  A(306) = RCT(306)*V(107)*V(164)
  A(307) = RCT(307)*V(107)*V(172)
  A(308) = RCT(308)*V(107)*V(165)
  A(309) = RCT(309)*V(39)
  A(310) = RCT(310)*V(87)*V(175)
  A(311) = RCT(311)*V(105)*V(164)
  A(312) = RCT(312)*V(105)*V(164)
  A(313) = RCT(313)*V(105)*V(172)
  A(314) = RCT(314)*V(86)*V(175)
  A(315) = RCT(315)*V(103)*V(164)
  A(316) = RCT(316)*V(103)*V(164)
  A(317) = RCT(317)*V(103)*V(172)
  A(318) = RCT(318)*V(86)*V(169)
  A(319) = RCT(319)*V(94)*V(175)
  A(320) = RCT(320)*V(133)*V(175)
  A(321) = RCT(321)*V(113)*V(173)
  A(322) = RCT(322)*V(113)*V(164)
  A(323) = RCT(323)*V(113)*V(165)
  A(324) = RCT(324)*V(113)*F(1)
  A(325) = RCT(325)*V(113)*F(1)
  A(326) = RCT(326)*V(111)*V(173)
  A(327) = RCT(327)*V(111)*V(164)
  A(328) = RCT(328)*V(111)*V(165)
  A(329) = RCT(329)*V(111)*F(1)
  A(330) = RCT(330)*V(111)*F(1)
  A(331) = RCT(331)*V(97)*V(173)
  A(332) = RCT(332)*V(97)*V(164)
  A(333) = RCT(333)*V(97)*V(165)
  A(334) = RCT(334)*V(97)*F(1)
  A(335) = RCT(335)*V(104)*V(173)
  A(336) = RCT(336)*V(104)*V(164)
  A(337) = RCT(337)*V(104)*V(165)
  A(338) = RCT(338)*V(104)*F(1)
  A(339) = RCT(339)*V(109)*V(173)
  A(340) = RCT(340)*V(109)*V(164)
  A(341) = RCT(341)*V(109)*V(165)
  A(342) = RCT(342)*V(109)*F(1)
  A(343) = RCT(343)*V(109)*F(1)
  A(344) = RCT(344)*V(101)*V(173)
  A(345) = RCT(345)*V(101)*V(164)
  A(346) = RCT(346)*V(101)*V(165)
  A(347) = RCT(347)*V(101)*F(1)
  A(348) = RCT(348)*V(101)*F(1)
  A(349) = RCT(349)*V(119)*V(173)
  A(350) = RCT(350)*V(119)*V(164)
  A(351) = RCT(351)*V(119)*V(165)
  A(352) = RCT(352)*V(119)*F(1)
  A(353) = RCT(353)*V(119)*F(1)
  A(354) = RCT(354)*V(85)*V(175)
  A(355) = RCT(355)*V(74)*V(175)
  A(356) = RCT(356)*V(79)*V(164)
  A(357) = RCT(357)*V(82)*V(164)
  A(358) = RCT(358)*V(79)*V(172)
  A(359) = RCT(359)*V(82)*V(172)
  A(360) = RCT(360)*V(79)
  A(361) = RCT(361)*V(79)
  A(362) = RCT(362)*V(159)
  A(363) = RCT(363)*V(67)
  A(364) = RCT(364)*V(163)
  A(365) = RCT(365)*V(108)
  A(366) = RCT(366)*V(108)
  A(367) = RCT(367)*V(168)
  A(368) = RCT(368)*V(138)
  A(369) = RCT(369)*V(125)
  A(370) = RCT(370)*V(125)
  A(371) = RCT(371)*V(125)
  A(372) = RCT(372)*V(151)
  A(373) = RCT(373)*V(160)
  A(374) = RCT(374)*V(160)
  A(375) = RCT(375)*V(160)
  A(376) = RCT(376)*V(156)
  A(377) = RCT(377)*V(155)
  A(378) = RCT(378)*V(98)
  A(379) = RCT(379)*V(71)
  A(380) = RCT(380)*V(63)
  A(381) = RCT(381)*V(59)
  A(382) = RCT(382)*V(60)
  A(383) = RCT(383)*V(75)
  A(384) = RCT(384)*V(70)
  A(385) = RCT(385)*V(68)
  A(386) = RCT(386)*V(88)
  A(387) = RCT(387)*V(81)
  A(388) = RCT(388)*V(89)
  A(389) = RCT(389)*V(84)
  A(390) = RCT(390)*V(90)
  A(391) = RCT(391)*V(83)
  A(392) = RCT(392)*V(80)
  A(393) = RCT(393)*V(92)
  A(394) = RCT(394)*V(174)
  A(395) = RCT(395)*V(49)
  A(396) = RCT(396)*V(131)
  A(397) = RCT(397)*V(118)
  A(398) = RCT(398)*V(112)
  A(399) = RCT(399)*V(139)
  A(400) = RCT(400)*V(133)
  A(401) = RCT(401)*V(78)
  A(402) = RCT(402)*V(110)
  A(403) = RCT(403)*V(94)
  A(404) = RCT(404)*V(37)
  A(405) = RCT(405)*V(74)
  A(406) = RCT(406)*V(150)
  A(407) = RCT(407)*V(44)*V(175)
  A(408) = RCT(408)*V(53)*V(175)
  A(409) = 1.7e-11*V(54)*V(175)
  A(410) = RCT(410)*V(127)*V(172)
  A(411) = RCT(411)*V(127)*V(164)
  A(412) = 2.3e-12*V(127)*V(166)
  A(413) = RCT(413)*V(127)*V(170)
  A(414) = RCT(414)*V(127)*V(167)
  A(415) = RCT(415)*V(77)*V(175)
  A(416) = 3.8e-12*V(77)*V(166)
  A(417) = 4.7e-11*V(76)*V(175)
  A(418) = 1.4e-11*V(76)*V(166)
  A(419) = 2e-11*V(93)*V(175)
  A(420) = 9.2e-18*V(93)*V(169)
  A(421) = 9.9e-11*V(93)*V(166)
  A(422) = RCT(422)*V(65)
  A(423) = RCT(423)*V(65)*V(175)
  A(424) = 2.4e-15*V(65)*V(166)
  A(425) = RCT(425)*V(124)*V(172)
  A(426) = RCT(426)*V(124)*V(164)
  A(427) = RCT(427)*V(124)*V(165)
  A(428) = 4.66e-12*V(51)*V(175)
  A(429) = RCT(429)*V(51)
  A(430) = RCT(430)*V(62)
  A(431) = 1.06e-12*V(62)*V(175)
  A(432) = 7e-12*V(129)*V(165)
  A(433) = RCT(433)*V(129)*V(164)
  A(434) = 2.3e-12*V(129)*V(166)
  A(435) = RCT(435)*V(129)*V(172)
  A(436) = 3.6e-12*V(66)*V(175)
  A(437) = RCT(437)*V(66)
  A(438) = 2.86e-13*V(102)*V(169)
  A(439) = 2.08e-12*V(102)*V(165)
  A(440) = RCT(440)*V(140)
  A(441) = 3.47e-12*V(140)*V(175)
  A(442) = 2.6e-12*V(140)*V(166)
  A(443) = RCT(443)*V(129)*V(170)
  A(444) = RCT(444)*V(124)*V(170)
  A(445) = 1.1e-10*V(95)*V(175)
  A(446) = 1.8e-15*V(95)*V(169)
  A(447) = 0.0017*V(95)
  A(448) = 1.1e-10*V(96)*V(175)
  A(449) = 1.8e-15*V(96)*V(169)
  A(450) = 0.0025*V(96)
  A(451) = RCT(451)*V(129)*V(171)
  A(452) = RCT(452)*V(124)*V(171)
  A(453) = RCT(453)*V(127)*V(171)
  A(454) = 9.9206e-11*F(2)
  A(455) = 1.1574e-11*F(2)
  A(456) = 3.9684e-12*F(2)
  A(457) = 3.3068e-13*F(2)
  A(458) = RCT(458)*F(2)
  A(459) = RCT(459)*F(2)
  A(460) = RCT(460)*F(2)
  A(461) = RCT(461)*F(2)
  A(462) = 1.1574e-05*V(128)
  A(463) = 1.1574e-05*V(108)
  A(464) = 1.1574e-05*V(116)
  A(465) = 1.1574e-05*V(159)
  A(466) = 1.1574e-05*V(45)
  A(467) = 1.1574e-05*V(142)
  A(468) = 1.1574e-05*V(78)
  A(469) = 1.1574e-05*V(137)
  A(470) = 1.1574e-05*V(44)
  A(471) = 1.1574e-05*V(1)
  A(472) = 1.1574e-05*V(47)
  A(473) = 1.1574e-05*V(36)
  A(474) = 1.1574e-05*V(149)
  A(475) = 1.1574e-05*V(113)
  A(476) = 1.1574e-05*V(101)
  A(477) = 1.1574e-05*V(173)
  A(478) = 1.1574e-05*V(2)
  A(479) = 1.1574e-05*V(3)
  A(480) = 1.1574e-05*V(37)
  A(481) = 1.1574e-05*V(87)
  A(482) = 1.1574e-05*V(79)
  A(483) = 1.1574e-05*V(105)
  A(484) = 1.1574e-05*V(72)
  A(485) = 1.1574e-05*V(94)
  A(486) = 1.1574e-05*V(147)
  A(487) = 1.1574e-05*V(63)
  A(488) = 1.1574e-05*V(119)
  A(489) = 1.1574e-05*V(138)
  A(490) = 1.1574e-05*V(125)
  A(491) = 1.1574e-05*V(73)
  A(492) = 1.1574e-05*V(155)
  A(493) = 1.1574e-05*V(85)
  A(494) = 1.1574e-05*V(153)
  A(495) = 1.1574e-05*V(143)
  A(496) = 1.1574e-05*V(91)
  A(497) = 1.1574e-05*V(50)
  A(498) = 1.1574e-05*V(117)
  A(499) = 1.1574e-05*V(58)
  A(500) = 1.1574e-05*V(172)
  A(501) = 1.1574e-05*V(74)
  A(502) = 1.1574e-05*V(82)
  A(503) = 1.1574e-05*V(84)
  A(504) = 1.1574e-05*V(48)
  A(505) = 1.1574e-05*V(46)
  A(506) = 1.1574e-05*V(41)
  A(507) = 1.1574e-05*V(121)
  A(508) = 1.1574e-05*V(4)
  A(509) = 1.1574e-05*V(55)
  A(510) = 1.1574e-05*V(157)
  A(511) = 1.1574e-05*V(98)
  A(512) = 1.1574e-05*V(110)
  A(513) = 1.1574e-05*V(150)
  A(514) = 1.1574e-05*V(5)
  A(515) = 1.1574e-05*V(135)
  A(516) = 1.1574e-05*V(123)
  A(517) = 1.1574e-05*V(120)
  A(518) = 1.1574e-05*V(90)
  A(519) = 1.1574e-05*V(114)
  A(520) = 1.1574e-05*V(112)
  A(521) = 1.1574e-05*V(134)
  A(522) = 1.1574e-05*V(139)
  A(523) = 1.1574e-05*V(130)
  A(524) = 1.1574e-05*V(148)
  A(525) = 1.1574e-05*V(6)
  A(526) = 1.1574e-05*V(156)
  A(527) = 1.1574e-05*V(131)
  A(528) = 1.1574e-05*V(107)
  A(529) = 1.1574e-05*V(111)
  A(530) = 1.1574e-05*V(132)
  A(531) = 1.1574e-05*V(92)
  A(532) = 1.1574e-05*V(136)
  A(533) = 1.1574e-05*V(49)
  A(534) = 1.1574e-05*V(7)
  A(535) = 1.1574e-05*V(8)
  A(536) = 1.1574e-05*V(167)
  A(537) = 1.1574e-05*V(168)
  A(538) = 1.1574e-05*V(109)
  A(539) = 1.1574e-05*V(151)
  A(540) = 1.1574e-05*V(104)
  A(541) = 1.1574e-05*V(86)
  A(542) = 1.1574e-05*V(103)
  A(543) = 1.1574e-05*V(170)
  A(544) = 1.1574e-05*V(57)
  A(545) = 1.1574e-05*V(52)
  A(546) = 1.1574e-05*V(152)
  A(547) = 1.1574e-05*V(80)
  A(548) = 1.1574e-05*V(160)
  A(549) = 1.1574e-05*V(118)
  A(550) = 1.1574e-05*V(97)
  A(551) = 1.1574e-05*V(38)
  A(552) = 1.1574e-05*V(146)
  A(553) = 1.1574e-05*V(164)
  A(554) = 1.1574e-05*V(165)
  A(555) = 1.1574e-05*V(166)
  A(556) = 1.1574e-05*V(100)
  A(557) = 1.1574e-05*V(34)
  A(558) = 1.1574e-05*V(169)
  A(559) = 1.1574e-05*V(175)
  A(560) = 1.1574e-05*V(67)
  A(561) = 1.1574e-05*V(39)
  A(562) = 1.1574e-05*V(9)
  A(563) = 1.1574e-05*V(144)
  A(564) = 1.1574e-05*V(70)
  A(565) = 1.1574e-05*V(40)
  A(566) = 1.1574e-05*V(145)
  A(567) = 1.1574e-05*V(115)
  A(568) = 1.1574e-05*V(71)
  A(569) = 1.1574e-05*V(133)
  A(570) = 1.1574e-05*V(10)
  A(571) = 1.1574e-05*V(162)
  A(572) = 1.1574e-05*V(174)
  A(573) = 1.1574e-05*V(158)
  A(574) = 1.1574e-05*V(75)
  A(575) = 1.1574e-05*V(59)
  A(576) = 1.1574e-05*V(60)
  A(577) = 1.1574e-05*V(163)
  A(578) = 1.1574e-05*V(161)
  A(579) = 1.1574e-05*V(154)
  A(580) = 1.1574e-05*V(88)
  A(581) = 1.1574e-05*V(81)
  A(582) = 1.1574e-05*V(89)
  A(583) = 1.1574e-05*V(106)
  A(584) = 1.1574e-05*V(68)
  A(585) = 1.1574e-05*V(53)
  A(586) = 1.1574e-05*V(141)
  A(587) = 1.1574e-05*V(83)
  A(588) = 1.1574e-05*V(54)
  A(589) = 1.1574e-05*V(42)
  A(590) = 1.1574e-05*V(11)
  A(591) = 1.1574e-05*V(64)
  A(592) = 1.1574e-05*V(12)
  A(593) = 1.1574e-05*V(122)
  A(594) = 1.1574e-05*V(69)
  A(595) = 1.1574e-05*V(13)
  A(596) = 1.1574e-05*V(14)
  A(597) = 1.1574e-05*V(15)
  A(598) = 1.1574e-05*V(16)
  A(599) = 1.1574e-05*V(35)
  A(600) = 1.1574e-05*V(43)
  A(601) = 1.1574e-05*V(126)
  A(602) = 1.1574e-05*V(171)
  A(603) = 1.1574e-05*V(61)
  A(604) = 1.1574e-05*V(56)
  A(605) = 1.1574e-05*V(99)
  A(606) = 1.1574e-05*V(76)
  A(607) = 1.1574e-05*V(77)
  A(608) = 1.1574e-05*V(65)
  A(609) = 1.1574e-05*V(129)
  A(610) = 1.1574e-05*V(66)
  A(611) = 1.1574e-05*V(62)
  A(612) = 1.1574e-05*V(102)
  A(613) = 1.1574e-05*V(140)
  A(614) = 1.1574e-05*V(124)
  A(615) = 1.1574e-05*V(93)
  A(616) = 1.1574e-05*V(96)
  A(617) = 1.1574e-05*V(17)
  A(618) = 1.1574e-05*V(127)

! Aggregate function
  Vdot(1) = -A(471)
  Vdot(2) = A(54)+A(74)+A(95)+A(100)+A(140)+A(149)+0.361*A(197)+0.5*A(208)+1.172*A(242)+0.858*A(243)+A(247)+2*A(250)+0.7&
              &*A(252)+2*A(253)+A(254)+A(255)+A(266)+A(272)+0.36*A(291)+0.63*A(292)+0.93*A(306)+0.93*A(307)+A(315)+0.5&
              &*A(317)+A(318)+A(319)+0.085*A(355)+0.4*A(402)+0.35*A(425)+A(426)+A(429)+A(431)+A(440)+A(444)+0.1*A(446)+0.1&
              &*A(449)+A(452)-A(478)
  Vdot(3) = 0.27*A(299)+0.21*A(301)-A(479)
  Vdot(4) = 0.781*A(194)-A(196)-A(508)
  Vdot(5) = 0.04*A(123)+0.04*A(124)+0.04*A(125)+0.04*A(126)+0.04*A(127)+0.04*A(128)+0.04*A(129)+0.04*A(130)+0.04*A(133)&
              &+0.04*A(134)+0.04*A(135)-A(514)
  Vdot(6) = 0.007*A(232)+0.007*A(233)+0.007*A(234)+0.007*A(235)+0.007*A(237)-A(238)+0.005*A(386)+0.005*A(387)+0.005&
              &*A(388)-A(525)
  Vdot(7) = -A(534)
  Vdot(8) = -A(535)
  Vdot(9) = -A(562)
  Vdot(10) = 0.35*A(304)-A(570)
  Vdot(11) = -A(590)
  Vdot(12) = -A(592)
  Vdot(13) = A(90)+0.25*A(126)+0.15*A(189)+A(221)+A(222)+A(358)+A(359)+A(360)+0.47*A(409)+0.45*A(445)+0.52*A(448)-A(595)
  Vdot(14) = -A(596)
  Vdot(15) = -A(597)
  Vdot(16) = -A(598)
  Vdot(17) = -A(617)
  Vdot(18) = A(204)
  Vdot(19) = A(216)
  Vdot(20) = A(407)
  Vdot(21) = A(408)
  Vdot(22) = A(409)
  Vdot(23) = A(200)+A(202)+A(369)+A(370)+A(371)+A(490)
  Vdot(24) = A(51)+A(56)+A(57)+A(474)
  Vdot(25) = A(95)+A(496)
  Vdot(26) = A(140)+A(464)
  Vdot(27) = A(201)+A(203)+A(372)+A(539)
  Vdot(28) = A(417)+A(418)+A(606)
  Vdot(29) = A(415)+A(416)+A(607)
  Vdot(30) = A(422)+A(423)+A(424)+A(608)
  Vdot(31) = A(436)+A(437)+A(610)
  Vdot(32) = A(440)+A(441)+A(442)+A(613)
  Vdot(33) = A(419)+A(420)+A(421)+A(615)
  Vdot(34) = -A(6)-A(13)+A(29)-A(557)
  Vdot(35) = -A(60)+A(454)-A(599)
  Vdot(36) = -A(103)-A(104)-A(473)
  Vdot(37) = A(361)-A(404)-A(480)
  Vdot(38) = A(12)-A(38)-A(40)-A(551)
  Vdot(39) = A(308)-A(309)-A(561)
  Vdot(40) = A(142)-A(143)-A(565)
  Vdot(41) = 0.5*A(237)-A(241)+0.1*A(297)-A(506)
  Vdot(42) = -A(42)+0.1*A(192)-A(589)
  Vdot(43) = A(71)+A(72)-A(74)+A(76)+A(77)+A(455)-A(600)
  Vdot(44) = -A(407)+A(459)-A(470)
  Vdot(45) = -A(107)-A(138)-A(466)
  Vdot(46) = 0.272*A(235)+0.32*A(236)-A(240)+0.033*A(298)-A(505)
  Vdot(47) = -A(101)-A(258)-A(472)
  Vdot(48) = 0.578*A(235)+0.68*A(236)-A(239)+0.067*A(298)-A(504)
  Vdot(49) = 0.41*A(188)-A(257)-A(395)-A(533)
  Vdot(50) = A(20)-A(27)-A(35)+A(440)-A(497)
  Vdot(51) = 0.65*A(425)-A(428)-A(429)
  Vdot(52) = A(47)-A(48)-A(58)-A(59)-A(545)
  Vdot(53) = -A(408)+A(460)-A(585)
  Vdot(54) = -A(409)+A(461)-A(588)
  Vdot(55) = -A(144)+A(146)+0.53*A(206)+0.5*A(218)+0.6*A(402)-A(509)
  Vdot(56) = A(67)-A(68)-A(78)-A(79)-A(604)
  Vdot(57) = A(44)-A(49)-A(50)-A(55)-A(544)
  Vdot(58) = A(24)-A(25)-A(39)-A(41)-A(499)
  Vdot(59) = A(186)-A(227)-A(381)-A(575)
  Vdot(60) = A(162)-A(228)-A(382)-A(576)
  Vdot(61) = A(62)-A(69)-A(70)-A(75)-A(603)
  Vdot(62) = A(427)-A(430)-A(431)-A(611)
  Vdot(63) = A(185)-A(226)-A(380)-A(487)
  Vdot(64) = -A(15)+A(57)+A(77)+0.09*A(174)+A(370)+0.19*A(389)-A(591)
  Vdot(65) = 0.06*A(408)+0.06*A(409)-A(422)-A(423)-A(424)-A(608)
  Vdot(66) = 0.15*A(425)+A(431)+A(435)-A(436)-A(437)-A(610)
  Vdot(67) = A(98)-A(99)+0.19*A(194)-A(363)-A(560)
  Vdot(68) = 0.75*A(126)+0.41*A(189)-A(230)-A(385)-A(584)
  Vdot(69) = -A(1)-A(2)-A(3)-A(4)-A(5)+A(6)+A(30)+A(32)+A(34)-A(594)
  Vdot(70) = A(187)-A(231)-A(384)-A(564)
  Vdot(71) = A(164)-A(225)-A(379)-A(568)
  Vdot(72) = 0.25*A(82)+0.25*A(166)-A(181)+A(184)+0.09*A(292)-A(484)
  Vdot(73) = -A(16)+A(19)-A(31)+A(324)+A(330)+A(334)+A(338)+A(343)+A(347)+A(352)-A(491)
  Vdot(74) = 0.5*A(114)-A(355)-A(405)-A(501)
  Vdot(75) = A(154)-A(229)+0.85*A(313)+0.5*A(317)-A(383)-A(574)
  Vdot(76) = 0.19*A(408)+0.15*A(409)-A(417)-A(418)-A(606)
  Vdot(77) = 0.54*A(407)-A(415)-A(416)-A(607)
  Vdot(78) = -A(93)-A(94)+0.85*A(156)+0.85*A(157)-A(401)-A(468)
  Vdot(79) = 0.5*A(114)-A(356)-A(358)-A(360)-A(361)-A(482)
  Vdot(80) = A(160)-A(246)-A(247)+A(329)-A(392)-A(547)
  Vdot(81) = 0.272*A(158)-A(233)-A(236)-A(387)-A(581)
  Vdot(82) = 0.365*A(355)-A(357)-A(359)-A(502)
  Vdot(83) = 0.1*A(159)-A(245)-A(391)-A(587)
  Vdot(84) = -A(244)+0.1*A(296)-A(389)-A(503)
  Vdot(85) = 0.44*A(242)+0.473*A(243)-A(354)-A(493)
  Vdot(86) = 0.17*A(267)+0.17*A(294)-A(314)-A(318)-A(541)
  Vdot(87) = 0.654*A(244)+0.29*A(267)+0.29*A(294)-A(310)-A(481)
  Vdot(88) = 0.628*A(158)-A(232)-A(235)-A(386)-A(580)
  Vdot(89) = 0.037*A(158)-A(234)-A(237)-A(388)-A(582)
  Vdot(90) = 0.5*A(300)+0.51*A(302)-A(303)-A(390)-A(518)
  Vdot(91) = -A(95)+0.134*A(197)+A(198)+0.5*A(208)+0.278*A(242)+0.142*A(243)+0.16*A(291)+0.057*A(299)+0.65*A(304)+0.07&
               &*A(306)+0.07*A(307)+A(325)+0.25*A(445)+0.5*A(446)+0.1*A(448)+0.2*A(449)-A(496)
  Vdot(92) = 0.85*A(161)+0.41*A(190)+0.5*A(218)-A(248)-A(249)-A(393)-A(531)
  Vdot(93) = 0.8*A(415)+0.727*A(417)-A(419)-A(420)-A(421)-A(615)
  Vdot(94) = 0.26*A(292)+0.21*A(299)+0.2*A(300)-A(319)-A(403)-A(485)
  Vdot(95) = 0.58*A(407)+0.7*A(408)+0.3*A(409)+0.18*A(415)+0.3*A(417)+0.3*A(418)+1.05*A(419)+A(420)+0.3*A(421)+A(440)&
               &+A(441)+A(442)-A(445)-A(446)-A(447)
  Vdot(96) = 0.35*A(408)+0.6*A(409)+0.2*A(417)+0.2*A(418)-A(448)-A(449)-A(450)-A(616)
  Vdot(97) = 0.044*A(213)-A(331)-A(332)-A(333)-A(334)-A(550)
  Vdot(98) = 0.96*A(129)+0.96*A(134)+A(163)-A(223)-A(224)-A(378)-A(511)
  Vdot(99) = A(63)+0.5*A(65)-A(73)+0.25*A(80)+0.25*A(81)+0.25*A(82)+0.25*A(83)+0.25*A(84)+0.5*A(85)+0.25*A(86)+0.25&
               &*A(87)+0.25*A(88)-A(605)
  Vdot(100) = A(145)-A(147)-A(193)-A(195)-A(556)
  Vdot(101) = 0.12*A(192)-A(344)-A(345)-A(346)-A(347)-A(348)-A(476)
  Vdot(102) = 0.06*A(415)+0.742*A(416)+0.073*A(417)+0.3*A(418)+0.3*A(419)+0.3*A(421)+A(432)+A(433)+A(434)+A(437)-A(438)&
                &-A(439)+A(443)+A(451)-A(612)
  Vdot(103) = A(314)-A(315)-A(316)-A(317)-A(542)
  Vdot(104) = 0.022*A(215)+0.055*A(259)-A(335)-A(336)-A(337)-A(338)-A(540)
  Vdot(105) = 0.03*A(113)+0.09*A(173)+0.113*A(266)-A(311)-A(312)-A(313)-A(483)
  Vdot(106) = 0.25*A(80)+0.25*A(81)+0.25*A(83)+0.25*A(84)+0.25*A(86)+0.25*A(87)+0.25*A(88)+0.25*A(152)+0.25*A(153)+0.25&
                &*A(168)+0.25*A(169)+0.25*A(171)+0.24*A(172)+0.25*A(174)+0.25*A(175)+0.15*A(176)+0.25*A(178)+0.25*A(180)&
                &-A(182)+0.3*A(252)-A(583)
  Vdot(107) = A(305)-A(306)-A(307)-A(308)+A(309)-A(528)
  Vdot(108) = 0.16*A(83)+0.75*A(87)+0.32*A(108)+A(119)-A(151)+0.16*A(168)+0.75*A(178)+0.791*A(228)+0.32*A(263)+A(270)&
                &+A(286)-A(365)-A(366)+A(382)+0.32*A(394)-A(463)
  Vdot(109) = 0.12*A(214)-A(339)-A(340)-A(341)-A(342)-A(343)-A(538)
  Vdot(110) = 0.96*A(125)+A(144)-A(146)-A(194)-A(402)-A(512)
  Vdot(111) = 0.044*A(213)-A(326)-A(327)-A(328)-A(329)-A(330)-A(529)
  Vdot(112) = 0.081*A(113)-A(291)-A(298)-A(398)-A(520)
  Vdot(113) = 0.12*A(192)+0.11*A(213)+0.12*A(214)+0.326*A(215)-A(321)-A(322)-A(323)-A(324)-A(325)-A(475)
  Vdot(114) = -A(204)-A(213)-A(216)-A(519)
  Vdot(115) = -A(191)-A(192)-A(199)+0.102*A(213)+A(373)-A(567)
  Vdot(116) = A(92)-A(140)+0.15*A(188)+0.5*A(208)+A(211)+A(212)+A(251)+A(262)+A(275)+A(276)+A(277)+A(278)+A(279)+A(280)&
                &+A(281)+A(282)+A(283)+A(284)+A(285)+A(286)+A(288)+0.08*A(306)+0.08*A(307)+A(348)+0.15*A(448)+0.3*A(449)&
                &-A(464)
  Vdot(117) = A(21)-A(28)-A(36)+A(52)+A(72)+A(97)+0.6*A(122)+A(138)+A(150)+A(167)+A(202)+A(203)+0.5*A(218)+A(256)+A(258)&
                &+0.41*A(291)+0.3*A(292)+A(295)+A(312)+A(316)+0.742*A(416)+0.5*A(418)+0.5*A(421)+A(424)+0.5*A(442)-A(498)
  Vdot(118) = 0.12*A(115)+0.21*A(291)+0.075*A(299)+0.06*A(300)+0.26*A(301)+0.16*A(302)-A(304)-A(397)-A(549)
  Vdot(119) = 0.055*A(259)-A(349)-A(350)-A(351)-A(352)-A(353)+A(420)+0.4*A(446)+0.25*A(449)-A(488)
  Vdot(120) = 0.655*A(132)-A(133)-A(134)-A(135)-A(517)
  Vdot(121) = A(239)+A(240)+A(241)-A(242)-A(243)-A(507)
  Vdot(122) = A(45)-A(53)+0.5*A(65)+0.25*A(130)+0.25*A(131)+0.25*A(135)+0.25*A(152)+0.25*A(153)+0.25*A(166)+0.25*A(168)&
                &+0.25*A(169)+0.5*A(170)+0.25*A(171)+0.24*A(172)+0.25*A(174)+0.25*A(175)+0.25*A(178)+0.25*A(179)+0.25*A(180)&
                &+0.24*A(292)+0.1*A(368)-A(593)
  Vdot(123) = 0.4*A(122)-A(127)-A(128)-A(129)-A(130)-A(516)
  Vdot(124) = A(423)+A(424)-A(425)-A(426)-A(427)+A(428)+A(430)-A(444)-A(452)-A(614)
  Vdot(125) = A(127)+0.94*A(128)+0.7*A(130)+0.7*A(131)+0.919*A(133)+0.7*A(135)+0.134*A(197)-A(200)-A(202)+0.085*A(242)&
                &+0.088*A(243)+0.216*A(267)+0.06*A(292)+0.216*A(294)+0.48*A(311)+0.072*A(313)+0.085*A(355)+A(356)+A(357)&
                &-A(369)-A(370)-A(371)+A(404)+0.25*A(406)+0.18*A(407)+0.12*A(408)+0.1*A(409)+0.06*A(415)+0.6*A(445)+A(446)&
                &-A(490)
  Vdot(126) = A(61)+A(63)+2*A(64)+0.5*A(65)+A(66)+A(70)-A(71)-A(72)+A(73)+A(75)-A(76)-A(77)+A(78)+0.75*A(80)+0.75*A(81)&
                &+0.75*A(82)+0.75*A(83)+0.75*A(84)+0.5*A(85)+0.75*A(86)+0.75*A(87)+0.75*A(88)+A(89)+A(90)+A(91)+A(92)+A(451)&
                &+A(452)+A(453)+A(457)-A(601)
  Vdot(127) = 0.46*A(407)+0.81*A(408)+0.85*A(409)-A(410)-A(411)-A(412)-A(413)-A(414)+0.14*A(415)+0.2*A(417)+0.2*A(418)&
                &+0.7*A(419)+0.2*A(421)-A(453)-A(618)
  Vdot(128) = -A(80)+0.03*A(83)+A(104)-A(105)+0.05*A(108)-A(152)+0.03*A(168)-A(186)-A(209)-A(211)+0.36*A(227)+0.05&
                &*A(263)+0.05*A(394)-A(462)
  Vdot(129) = A(422)+0.2*A(425)+A(426)+A(429)-A(432)-A(433)-A(434)-A(435)+A(436)+A(438)-A(443)+A(444)-A(451)+A(452)&
                &-A(609)
  Vdot(130) = 0.9*A(297)-0.943*A(299)-A(300)+0.193*A(303)-A(523)
  Vdot(131) = 0.15*A(116)+0.32*A(291)+0.019*A(299)+0.01*A(300)+0.44*A(301)+0.28*A(302)-A(305)-A(396)-A(527)
  Vdot(132) = -A(118)-A(161)-A(177)+A(217)-A(272)-A(283)-A(530)
  Vdot(133) = 0.5*A(118)+0.96*A(123)+0.96*A(124)+0.94*A(127)+0.9*A(128)+0.66*A(130)+0.3*A(131)+0.894*A(133)+0.66*A(135)&
                &+0.075*A(161)+0.375*A(177)+0.791*A(225)+0.5*A(272)+0.42*A(291)+0.38*A(292)+0.3*A(299)+0.23*A(300)-A(320)&
                &-A(400)-A(569)
  Vdot(134) = 0.9*A(298)-A(301)-A(302)+0.193*A(303)-A(521)
  Vdot(135) = 0.6*A(122)-A(123)-A(124)-A(125)-A(126)+0.345*A(132)-A(515)
  Vdot(136) = A(249)-A(250)-A(251)-A(252)-2*A(253)-A(254)-A(255)-A(256)-A(532)
  Vdot(137) = 0.09*A(83)-A(87)+A(103)+0.18*A(108)-A(119)-A(162)+0.09*A(168)-A(178)+0.209*A(228)+0.18*A(263)-A(270)&
                &-A(286)+0.18*A(394)-A(469)
  Vdot(138) = 0.53*A(115)+0.578*A(159)+0.13*A(174)+0.36*A(175)-A(197)-A(198)+0.025*A(242)+0.088*A(243)+0.5*A(259)+0.234&
                &*A(267)+0.6*A(268)+0.24*A(292)+0.234*A(294)+0.9*A(296)+0.3*A(299)+0.23*A(300)+0.09*A(301)+0.06*A(302)+0.52&
                &*A(311)+0.078*A(313)+A(349)+A(350)+A(351)+A(352)-A(368)+0.26*A(389)+0.7*A(391)+A(397)+0.11*A(405)+0.5&
                &*A(406)-A(489)
  Vdot(139) = 0.009*A(113)+0.5*A(136)+0.25*A(179)-A(292)-A(297)-A(399)-A(522)
  Vdot(140) = 0.258*A(416)+0.5*A(418)+0.5*A(421)+A(439)-A(440)-A(441)-A(442)-A(613)
  Vdot(141) = -A(115)-A(159)-A(175)+A(205)+0.054*A(213)+0.209*A(245)-A(268)-A(280)-A(586)
  Vdot(142) = -A(85)+A(93)-A(111)+A(151)-A(156)-A(170)-A(264)-A(276)-A(467)
  Vdot(143) = -A(174)+0.143*A(235)+0.493*A(237)+0.346*A(244)-A(267)-A(279)+A(293)-A(294)-A(295)-A(296)-A(495)
  Vdot(144) = -A(81)-A(106)-A(153)-A(187)+A(191)-A(210)-A(212)+0.209*A(231)-A(563)
  Vdot(145) = -A(88)-A(137)-A(164)-A(180)+A(199)+0.209*A(225)-A(274)-A(285)-A(566)
  Vdot(146) = -A(145)+A(147)-A(149)-A(190)-A(220)-A(222)+A(248)-A(290)-A(552)
  Vdot(147) = -A(82)+0.16*A(83)+A(89)+A(101)-A(102)+0.32*A(108)+A(148)-A(166)+0.16*A(168)-2*A(183)-2*A(184)-A(185)+0.44&
                &*A(189)+A(219)+0.36*A(226)+A(258)+0.32*A(263)-A(287)-A(288)+A(289)+A(364)+0.85*A(367)+0.32*A(394)-A(486)
  Vdot(148) = -A(86)-A(112)-A(157)+A(165)+A(167)-A(171)-A(265)-A(277)-A(524)
  Vdot(149) = A(43)+A(45)+2*A(46)+A(50)-A(51)-A(52)+A(53)+A(55)-A(56)-A(57)+A(58)+0.5*A(65)+A(66)+0.5*A(81)+0.2*A(84)&
                &+0.5*A(88)+0.05*A(96)+A(106)+0.57*A(110)+A(111)+0.82*A(113)+0.35*A(115)+0.122*A(116)+0.5*A(118)+0.07*A(120)&
                &+0.07*A(121)+0.75*A(130)+0.75*A(131)+0.75*A(135)+0.07*A(136)+A(137)+0.5*A(149)+0.75*A(152)+1.25*A(153)+0.15&
                &*A(156)+0.063*A(158)+0.102*A(159)+0.075*A(161)+0.75*A(166)+0.75*A(168)+0.95*A(169)+0.5*A(170)+0.75*A(171)&
                &+1.22*A(172)+0.75*A(173)+0.95*A(174)+0.89*A(175)+1.255*A(176)+0.625*A(177)+0.75*A(178)+0.785*A(179)+1.25&
                &*A(180)+0.59*A(190)+0.5*A(192)+0.19*A(194)+A(195)+0.732*A(197)+A(210)+0.456*A(213)+0.6*A(214)+0.12*A(215)&
                &+A(219)+2*A(220)+A(221)+A(222)+0.375*A(242)+0.375*A(243)+A(252)+0.445*A(259)+A(261)+A(262)+A(264)+0.747&
                &*A(266)+0.4*A(268)+0.143*A(269)+0.39*A(271)+0.5*A(272)+0.07*A(273)+A(274)+A(290)+0.62*A(291)+0.2*A(292)&
                &+0.15*A(299)+0.07*A(300)+0.7*A(301)+0.44*A(302)+0.35*A(304)+0.08*A(306)+0.08*A(307)+A(319)+A(321)+A(322)&
                &+A(323)+A(324)+A(354)+0.55*A(355)+0.9*A(368)+A(371)+A(374)+A(376)+A(377)+A(384)+0.71*A(386)+0.71*A(387)&
                &+0.71*A(388)+0.3*A(391)+0.5*A(392)+A(393)+0.5*A(396)+A(400)+A(401)+0.4*A(402)+A(403)+A(406)+A(413)+A(443)&
                &+A(444)-A(474)
  Vdot(150) = 0.7*A(120)+0.7*A(121)-A(122)-A(131)-A(132)+1.2*A(136)+0.6*A(179)+0.7*A(273)-A(406)-A(513)
  Vdot(151) = 0.5*A(85)+A(94)+0.35*A(115)+0.122*A(116)+0.5*A(118)+0.102*A(159)+0.075*A(161)+0.5*A(170)+0.29*A(174)+0.14&
                &*A(175)+0.255*A(176)+0.375*A(177)-A(201)-A(203)+A(207)+0.5*A(214)+0.88*A(215)+0.085*A(242)+0.088*A(243)+0.5&
                &*A(259)+0.234*A(267)+0.4*A(268)+0.143*A(269)+0.5*A(272)+A(276)+0.234*A(294)+0.9*A(296)+0.65*A(304)+0.07&
                &*A(306)+0.07*A(307)+0.52*A(311)+0.078*A(313)+A(320)+A(335)+A(336)+A(337)+A(338)+A(339)+A(340)+A(341)+A(343)&
                &+0.5*A(354)+0.55*A(355)+A(356)+A(357)-A(372)+0.58*A(389)+0.3*A(391)+A(396)+A(404)+0.5*A(405)+0.25*A(406)&
                &+0.12*A(408)+0.2*A(409)+0.6*A(448)+A(449)-A(539)
  Vdot(152) = -A(116)-A(117)-A(160)-A(176)+0.47*A(206)+0.5*A(218)+A(246)-A(269)-A(281)-A(546)
  Vdot(153) = 0.058*A(113)+0.23*A(120)+0.23*A(121)+0.23*A(136)+0.3*A(172)+1.11*A(173)+0.115*A(179)+0.243*A(232)+0.513&
                &*A(233)+0.743*A(234)-A(259)+0.14*A(266)+0.23*A(273)-A(293)+0.275*A(386)+0.275*A(387)+0.275*A(388)+A(398)&
                &+A(399)-A(494)
  Vdot(154) = -A(113)-A(114)-A(158)-A(172)-2*A(173)+A(204)+0.75*A(232)+0.48*A(233)+0.25*A(234)-A(266)-A(278)-A(579)
  Vdot(155) = 0.16*A(81)+0.2*A(85)+0.728*A(116)+A(117)+0.16*A(153)+0.2*A(170)+0.18*A(174)+0.595*A(176)+A(193)+0.03&
                &*A(194)-A(207)-A(208)+0.65*A(212)+0.791*A(231)+0.085*A(242)+0.117*A(243)+A(247)+A(250)+0.7*A(252)+2*A(253)&
                &+A(254)+A(255)+0.216*A(267)+0.857*A(269)+0.42*A(292)+0.216*A(294)+0.27*A(299)+0.2*A(300)+0.09*A(301)+0.06&
                &*A(302)+0.85*A(306)+0.85*A(307)+0.5*A(310)+0.48*A(311)+0.072*A(313)+A(353)-A(377)+0.36*A(389)+A(392)+A(396)&
                &+0.39*A(405)+0.5*A(406)-A(492)
  Vdot(156) = 0.344*A(113)+0.035*A(120)+0.035*A(121)+0.035*A(136)+0.025*A(158)+0.18*A(172)+0.29*A(173)+0.018*A(179)&
                &-A(206)+0.3*A(213)-A(215)-A(217)-A(218)+0.294*A(266)+0.035*A(273)+A(326)+A(327)+A(328)+A(330)-A(376)+0.285&
                &*A(386)+0.285*A(387)+0.285*A(388)-A(526)
  Vdot(157) = -A(120)-A(121)-2*A(136)-A(163)-A(179)+A(216)+0.36*A(224)-A(273)-A(284)-A(510)
  Vdot(158) = -A(83)+0.15*A(84)+A(107)-A(108)-A(109)+A(138)-A(154)-A(168)+0.15*A(169)+0.209*A(229)-A(263)+0.3*A(271)&
                &-A(275)+0.25*A(450)-A(573)
  Vdot(159) = 0.5*A(81)+0.75*A(82)+0.16*A(83)+0.38*A(84)+0.5*A(86)+0.5*A(88)-A(96)-A(97)+A(102)+A(106)+0.32*A(108)+0.86&
                &*A(110)+0.93*A(112)+A(137)+0.5*A(153)+0.15*A(157)+0.75*A(166)+0.16*A(168)+0.38*A(169)+0.5*A(171)+0.5*A(180)&
                &+A(181)+2*A(183)+A(184)+0.5*A(192)+A(210)+0.1*A(214)+0.64*A(226)+0.32*A(263)+A(265)+0.75*A(271)+A(274)&
                &+A(287)+A(288)+A(344)+A(345)+A(346)+A(347)-A(362)+A(380)+A(384)+A(385)+0.32*A(394)-A(465)
  Vdot(160) = 0.476*A(113)+0.035*A(120)+0.035*A(121)+0.035*A(136)+0.038*A(158)+0.28*A(172)+0.45*A(173)+0.018*A(179)&
                &-A(205)+0.2*A(213)-A(214)+0.453*A(266)+0.035*A(273)+A(331)+A(332)+A(333)+A(334)-A(373)-A(374)-A(375)+0.425&
                &*A(386)+0.425*A(387)+0.425*A(388)+0.15*A(450)-A(548)
  Vdot(161) = -A(89)-A(90)+A(141)-A(142)+A(143)-A(148)+A(150)-A(189)-A(219)-A(221)+A(230)-A(289)+0.15*A(367)+A(375)+0.33&
                &*A(445)+0.25*A(447)+0.33*A(448)-A(578)
  Vdot(162) = -A(84)-A(110)+A(139)-A(155)-A(169)-A(271)-A(282)+0.612*A(303)+0.5*A(441)+0.5*A(442)-A(571)
  Vdot(163) = 0.75*A(80)+0.09*A(81)+0.07*A(83)+0.54*A(84)+0.25*A(88)+A(105)+0.14*A(108)+0.57*A(110)-A(141)-A(150)+0.75&
                &*A(152)+0.09*A(153)+0.033*A(159)+0.07*A(168)+0.54*A(169)+0.25*A(180)+A(182)+A(209)+A(211)+0.35*A(212)+0.64&
                &*A(227)+0.791*A(229)+0.13*A(263)+0.09*A(267)+0.57*A(271)+A(282)+A(283)+A(284)+A(285)+0.09*A(294)+A(315)+0.5&
                &*A(317)-A(364)+A(378)+A(379)+A(381)+A(383)+A(390)+0.13*A(394)+A(447)+A(450)-A(577)
  Vdot(164) = -A(3)+A(4)-A(7)-2*A(9)-A(10)+A(11)-A(20)-A(23)+A(32)+A(33)+A(35)-A(43)-A(61)-A(100)-A(102)-A(105)-A(106)&
                &-A(108)-A(109)-A(110)-A(111)-A(112)-A(113)-A(115)-A(116)-A(118)-A(119)-A(120)-A(124)-A(128)-A(133)-A(137)&
                &-A(148)-A(149)-A(243)-A(255)-A(256)-A(294)-A(295)-A(299)-A(301)-A(306)-A(311)-A(312)-A(315)-A(316)-A(322)&
                &-A(327)-A(332)-A(336)-A(340)-A(345)-A(350)-A(356)-A(357)-A(411)-A(426)-A(433)+A(458)-A(553)
  Vdot(165) = A(3)-A(4)-A(5)+A(7)-A(8)+2*A(9)+2*A(10)-A(12)-A(21)+A(22)+A(23)-A(24)+A(25)+A(26)+A(27)-A(32)+A(34)+A(36)&
                &+2*A(37)+A(38)+0.95*A(39)+A(40)+A(41)+A(43)-A(47)+A(48)+A(59)+A(61)-A(67)+A(68)+A(79)+A(84)+A(88)-A(98)&
                &+A(99)+A(100)+A(102)+A(105)+A(106)+A(108)+2*A(110)+A(111)+0.93*A(112)+0.91*A(113)+0.88*A(115)+0.85*A(116)&
                &+1.5*A(118)+A(119)+1.3*A(120)+1.3*A(121)+A(123)+A(124)-A(125)+0.96*A(126)+2*A(127)+1.88*A(128)+0.7*A(130)&
                &+0.7*A(131)+0.3*A(136)+2*A(137)-A(142)+A(143)-A(144)-A(145)+A(146)+A(147)+A(148)+A(149)+0.925*A(161)+A(169)&
                &+0.375*A(177)+0.15*A(179)+A(180)+A(193)+A(223)+A(243)+A(255)+A(271)+0.5*A(272)+0.3*A(273)+A(274)+A(282)&
                &+A(283)+A(284)+A(285)+0.06*A(292)+A(294)+0.1*A(297)+0.1*A(298)+0.79*A(299)+0.88*A(301)+0.06*A(302)+1.85&
                &*A(306)+0.85*A(307)-A(308)+A(309)+A(311)+A(315)+A(319)+A(320)+A(322)-A(323)+A(327)-A(328)+A(332)-A(333)&
                &+A(336)-A(337)+A(340)-A(341)+A(345)-A(346)+A(350)-A(351)+A(356)+A(357)+0.7*A(363)+A(378)+A(379)+A(390)&
                &+A(394)+A(396)+A(397)+A(398)+A(399)+A(400)+0.6*A(402)+A(403)+A(406)+A(411)+A(412)+A(426)-A(427)+A(430)&
                &+A(431)-A(432)+A(433)+A(434)-A(439)+0.5*A(441)+A(442)-A(554)
  Vdot(166) = A(5)+A(8)-A(10)-A(11)-A(12)-A(22)-A(26)+A(28)-A(33)-A(34)-2*A(37)+A(38)+0.05*A(39)+A(40)-A(52)+A(58)-A(72)&
                &+A(78)-A(97)-A(121)-A(122)-A(123)-A(127)-A(138)-A(150)-A(167)+0.81*A(194)+A(195)-A(199)-A(202)-A(203)&
                &-A(216)-A(217)-A(218)-A(258)+A(304)+0.15*A(306)+0.15*A(307)+A(323)+A(328)+A(333)+A(337)+A(341)+A(346)&
                &+A(351)+0.3*A(363)+0.4*A(402)-A(412)-A(416)-A(418)-A(421)-A(424)+A(432)-A(434)-A(442)-A(555)
  Vdot(167) = 0.3*A(85)+0.5*A(86)-A(91)-A(92)+0.95*A(96)+A(97)-A(98)+A(99)-A(100)+A(111)+0.93*A(112)+0.53*A(115)+0.5&
                &*A(149)+0.15*A(156)+0.15*A(157)+0.578*A(159)+0.3*A(170)+0.5*A(171)+0.36*A(175)-A(188)+A(195)+A(201)+A(203)&
                &-A(209)-A(210)-A(211)-A(212)+0.105*A(213)+0.28*A(214)+0.098*A(215)+A(220)-A(250)-A(251)+A(257)+0.445*A(259)&
                &-2*A(260)-A(261)-A(262)-A(263)-A(266)-A(267)-0.4*A(268)-A(269)-A(270)-A(271)-A(272)-A(273)-A(274)-A(275)&
                &-A(276)-A(277)-A(278)-A(279)-A(280)-A(281)-A(282)-A(283)-A(284)-A(285)-A(286)-A(287)-A(288)-A(289)+0.1&
                &*A(292)+A(342)+0.5*A(354)+0.085*A(355)+0.12*A(362)+0.7*A(363)+A(365)+0.85*A(367)+A(372)+A(374)+A(376)&
                &+A(377)+0.7*A(391)+A(393)+A(397)+A(401)+0.4*A(402)+0.11*A(405)-A(414)-A(536)
  Vdot(168) = 0.35*A(83)+0.25*A(86)+0.19*A(108)+0.187*A(159)-A(165)-A(167)+0.35*A(168)+0.25*A(171)+0.25*A(174)+0.25&
                &*A(175)+A(223)+0.791*A(245)+A(251)+0.19*A(263)+A(275)+A(277)+A(278)+A(279)+A(280)+A(281)+0.5*A(310)+A(318)&
                &-A(367)+0.19*A(394)-A(537)
  Vdot(169) = A(1)-A(2)-A(7)-A(8)-A(14)-A(17)-A(29)-A(30)+0.25*A(126)-A(131)+0.15*A(188)+0.15*A(189)+0.15*A(190)-A(192)&
                &-A(195)-A(213)-A(214)-A(215)-A(259)-A(291)-A(292)-A(318)-A(420)+0.15*A(425)-A(438)-A(446)-A(449)+A(456)&
                &-A(558)
  Vdot(170) = A(42)-A(43)-A(44)-2*A(45)-2*A(46)-A(47)+A(48)+A(49)+A(59)-A(65)-A(66)+0.09*A(83)+A(91)+A(100)+0.19*A(108)&
                &-A(130)-A(135)+A(140)+0.5*A(149)-A(152)-A(153)-A(166)-0.91*A(168)-A(169)-A(170)-A(171)-A(172)-A(174)-A(175)&
                &-A(176)-A(177)-A(178)-A(179)-A(180)+0.44*A(188)+0.39*A(190)+0.28*A(192)+0.5*A(208)+A(209)+A(210)+0.051&
                &*A(213)-A(219)-A(220)-A(221)-A(222)+A(250)-A(252)+2*A(260)-A(262)+A(263)+A(264)+A(265)+A(266)+A(267)+A(268)&
                &+A(269)+A(270)+A(271)+A(272)+A(273)+A(274)+A(287)+A(289)+A(290)+A(342)+0.88*A(362)+0.3*A(363)+A(365)+2&
                &*A(366)+0.15*A(367)+A(375)+0.18*A(394)+A(395)+A(400)-A(413)+A(414)-A(443)-A(444)-A(543)
  Vdot(171) = A(60)-A(61)-A(62)-2*A(63)-2*A(64)-A(65)-A(66)-A(67)+A(68)+A(69)+A(79)-A(80)-A(81)-A(82)-A(83)-A(84)-A(85)&
                &-A(86)-A(87)-A(88)-A(89)-A(90)-A(91)-A(92)-A(451)-A(452)-A(453)-A(602)
  Vdot(172) = A(14)+A(15)+A(16)-A(17)-A(18)-2*A(19)+A(22)-A(23)-A(24)-A(26)+0.95*A(39)+A(41)+A(43)-A(44)+2*A(46)+A(51)&
                &+A(52)+A(53)+A(54)+A(55)+2*A(56)+A(58)+A(61)-A(62)+2*A(64)+2*A(66)+A(71)+A(72)+A(73)+A(74)+A(75)+2*A(76)&
                &+A(78)+A(80)+A(81)+A(82)+0.64*A(83)+0.5*A(84)+0.3*A(85)+0.5*A(86)+A(87)+0.5*A(88)+A(89)+A(91)+A(95)+0.05&
                &*A(96)+A(102)+A(105)+A(106)+0.27*A(108)+0.82*A(113)+0.35*A(115)+0.85*A(116)+A(119)+0.8*A(120)+0.8*A(121)&
                &+A(123)+A(124)-A(126)-A(129)+0.5*A(130)+0.5*A(131)+0.934*A(133)-A(134)+1.2*A(135)+0.8*A(136)+A(152)+A(153)&
                &-A(154)-A(155)-A(156)-A(157)-0.937*A(158)-0.898*A(159)-A(160)-0.925*A(161)-A(162)-A(163)-A(164)+A(166)+0.64&
                &*A(168)+0.5*A(169)+0.3*A(170)+0.5*A(171)+1.1*A(172)+0.91*A(173)+A(174)+0.64*A(175)+1.7*A(176)+0.375*A(177)&
                &+A(178)+0.9*A(179)+0.5*A(180)+A(181)+A(182)+2*A(183)-A(185)-A(186)-A(187)-A(188)-A(189)-A(190)+0.28*A(192)&
                &+0.19*A(194)+0.773*A(197)+A(200)+A(202)+A(207)+A(209)+A(210)+0.237*A(213)+0.28*A(214)+0.141*A(215)+0.5&
                &*A(218)+A(219)+A(220)-0.175*A(242)+0.825*A(243)+0.7*A(252)-A(254)+0.445*A(259)+A(261)+0.27*A(263)+0.887&
                &*A(266)+A(267)+0.4*A(268)+0.85*A(269)+A(270)+0.5*A(272)+0.8*A(273)+A(287)+0.05*A(291)+A(294)-0.1*A(296)&
                &+0.61*A(299)-0.5*A(300)+0.69*A(301)-0.56*A(302)-A(307)+A(310)+A(311)-0.85*A(313)+A(315)-0.5*A(317)+A(318)&
                &+0.5*A(354)-A(358)-A(359)+A(362)+A(364)+1.73*A(368)+2*A(369)+A(372)+A(374)+A(376)+A(377)+A(378)+A(379)&
                &+A(380)+A(381)+A(382)+A(383)+A(384)+A(385)+0.985*A(386)+0.985*A(387)+0.985*A(388)+A(389)+A(390)+0.3*A(391)&
                &+A(392)+0.27*A(394)+A(396)+A(398)+A(399)+A(403)+0.89*A(405)+A(406)+0.54*A(407)+0.19*A(408)+0.15*A(409)&
                &+A(411)+A(412)+2*A(413)+A(414)+0.8*A(415)+0.727*A(417)+A(420)+A(422)-A(425)-A(435)+A(440)+A(443)+A(444)&
                &+0.33*A(445)+0.1*A(446)+0.25*A(447)+0.33*A(448)+0.1*A(449)+0.4*A(450)+A(451)+A(452)+2*A(453)-A(500)
  Vdot(173) = A(51)+A(52)-A(54)+A(56)+A(57)+0.05*A(96)+0.728*A(116)+A(117)+0.5*A(118)+A(123)+A(124)+0.5*A(149)+0.075&
                &*A(161)+0.29*A(173)+0.33*A(174)+0.595*A(176)+0.375*A(177)+0.39*A(190)+0.56*A(192)+A(193)+0.03*A(194)+0.505&
                &*A(197)+A(198)+2*A(200)+A(201)+2*A(202)+A(203)+0.5*A(208)+0.495*A(213)+0.56*A(214)+0.625*A(215)+0.5*A(218)&
                &+0.6*A(242)+0.678*A(243)+0.89*A(259)+0.09*A(267)+0.857*A(269)+0.5*A(272)+0.21*A(291)+0.09*A(294)+1.5*A(310)&
                &-A(321)-A(326)-A(331)-A(335)-A(339)-A(344)-A(349)+0.5*A(354)+0.55*A(355)+A(360)+0.88*A(362)+A(364)+A(366)&
                &+A(368)+2*A(369)+2*A(370)+A(371)+A(372)+A(373)+A(374)+A(376)+0.67*A(389)+0.5*A(392)+0.5*A(396)+A(400)&
                &+A(403)+1.89*A(405)+0.08*A(407)+0.1*A(408)+0.12*A(409)+A(422)+A(440)+0.25*A(445)+0.6*A(446)+0.25*A(447)&
                &+0.25*A(448)+0.6*A(449)+0.4*A(450)-A(477)
  Vdot(174) = A(109)+0.07*A(112)-A(139)+A(155)+0.25*A(177)+0.64*A(224)-A(394)-A(572)
  Vdot(175) = 2*A(13)-A(14)-A(15)-A(16)+A(17)-A(18)-A(20)-A(21)-A(22)+A(23)-A(25)+A(26)-A(27)-A(28)+2*A(31)+A(35)+A(36)&
                &+0.05*A(39)-A(42)-A(49)-A(51)-A(53)-A(54)+A(55)-A(60)-A(69)-A(71)-A(73)-A(74)+A(75)-A(93)-A(95)-A(96)&
                &-A(101)-A(103)-A(104)-A(107)+A(117)+0.5*A(118)-A(132)-A(139)-A(140)-A(141)-A(151)+0.15*A(156)+0.15*A(157)&
                &+0.063*A(158)+0.68*A(159)+0.15*A(161)-A(165)-A(181)-A(182)+0.44*A(188)+0.44*A(189)+0.44*A(190)-A(191)+0.36&
                &*A(192)-A(193)-A(194)-0.773*A(197)-A(200)-A(201)-A(204)-A(205)-A(206)-A(207)+0.243*A(213)+0.36*A(214)+0.239&
                &*A(215)-0.36*A(224)-0.209*A(225)-0.36*A(226)-0.36*A(227)-0.209*A(228)-0.209*A(229)-A(230)-0.209*A(231)&
                &-0.875*A(232)-0.74*A(233)-0.625*A(234)-0.15*A(235)-0.5*A(237)-A(239)-A(240)-A(241)+1.125*A(242)+0.125&
                &*A(243)-0.346*A(244)-0.209*A(245)-A(246)-A(248)-A(249)+A(250)+0.7*A(252)+2*A(253)+2*A(254)+A(255)-A(257)&
                &+0.89*A(259)+0.05*A(291)+0.37*A(292)-A(293)+0.9*A(296)-A(297)-A(298)+0.5*A(300)+0.5*A(302)-0.388*A(303)&
                &-A(304)-A(305)+A(307)-A(310)+0.15*A(313)-A(314)+0.5*A(317)+A(318)-A(319)-A(320)-A(354)-0.365*A(355)+A(356)&
                &+A(357)+A(360)+0.07*A(368)+A(378)+A(379)+A(380)+A(381)+A(382)+A(383)+A(384)+A(385)+0.985*A(386)+0.985&
                &*A(387)+0.985*A(388)+A(389)+A(390)+A(391)+A(392)+A(393)+A(395)+A(401)+2*A(404)+2*A(405)-A(407)-A(408)&
                &-A(409)+A(410)-A(415)-A(417)-A(419)+A(420)-A(423)+0.2*A(425)-A(428)+A(429)-A(431)-A(436)-A(441)-0.75*A(445)&
                &+0.1*A(446)-0.75*A(448)+0.1*A(449)-A(559)
      
END SUBROUTINE Fun

! End of Fun function
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~



END MODULE aromatics_kpp_Function

